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Abstract. We present the results of a near-infrared adaptive optics survey with the aim to detect close companions to Hipparcos 
members in the three subgroups of the nearby OB association Sco OB2: Upper Scorpius (US), Upper Centaurus Lupus (UCL) 
and Lower Centaurus Crux (LCC). We have targeted 199 A-type and late B-type stars in the K s band, and a subset also in 
the J and H band. We find 151 stellar components other than the target stars. A brightness criterion is used to separate these 
components into 77 background stars and 74 candidate physical companion stars. Out of these 74 candidate companions, 41 
have not been reported before (14 in US; 13 in UCL; 14 in LCC). The angular separation between primaries and observed 
companion stars ranges from 0.22" to 12.4". At the mean distance of Sco OB2 (130 pc) this corresponds to a projected 
separation of 28.6 AU to 1612 AU. Absolute magnitudes are derived for all primaries and observed companions using the 
parallax and interstellar extinction for each star individually. For each object we derive the mass from K s , assuming an age 
of 5 Myr for the US subgroup, and 20 Myr for the UCL and LCC subgroups. Companion star masses range from 0.10 M G to 
3.0 M Q . The mass ratio distribution follows f(q) = q~ r with T = 0.33, which excludes random pairing. No close (p < 3.75") 
companion stars or background stars are found in the magnitude range 12 mag < Ks < 14 mag. The lack of stars with these 
properties cannot be explained by low-number statistics, and may imply a lower limit on the companion mass of ~ 0.1 M Q . 
Close stellar components with K s > 14 mag are observed. If these components are very low-mass companion stars, a gap in the 
companion mass distribution might be present. The small number of close low-mass companion stars could support the embryo- 
ejection formation scenario for brown dwarfs. Our findings are compared with and complementary to visual, spectroscopic, and 
astrometric data on binarity in Sco OB2. We find an overall companion star fraction of 0.52 in this association. This is a lower 
limit since the data from the observations and from literature are hampered by observational biases and selection effects. This 
paper is the first step toward our goal to derive the primordial binary population in Sco OB2. 

Key words, binaries: visual - binaries: general - stars: formation - associations - individual: Sco OB2 



1. Introduction 

Duplicity and multiplicity properties of newly born stars are 
among the most importan t clues to understanding the process 
of star formation (Blaauwl ll99ll) . Observations of star forming 
regions over the past two decades have revealed two important 
facts: (1) prac tically all (70-90%) stars form in clusters (e.g., 
lLada & LadaL l2003h a nd, (2) w i thin th ese clusters most stars 
are formed in binaries (Mathieu, 1994). Consequently, the star 
formation community has shifted its attention toward under- 
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standing the formation of multiple systems — from binaries to 
star clusters — by means of both observations and theory. 

The observational progress in studies of very young em- 
bedded as well as expo sed star clusters is e xtensively summa- 
rized in the review bv lLada&Ladal §003). On the theoreti- 
cal side the numerical simulations of cluster formation have 
become increasingly sophisticated, covering the very earliest 
phases of the develop ment of massive dense cores in gian t 
molecular clouds (e.g., Klessen, Heitsch, & M ac Lowt |2000), 
the subsequent clustered formation of stars and binaries (e.g., 
iBate. Bonnell. & Brommll2003l) . as well as the early evolution 
of the binary population during the phase of gas expulsion from 
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the newly formed cluster dKroupa. Aarseth. & Hurlevl |2001 ). 
At the same time numerical simulations of older exposed clus- 
ters have become more realistic by incorporating detailed stel- 
lar and binary evolution eff ects in N-body simulations (e.g., 
IPortegies Zwart et all 1200 ll) . This has led to the creation of a 
number of research networks that aim at synthesizing the mod- 
eling and observing efforts into a single framework which cov- 
ers all the stages from the formation of a star cluster to its even- 
tual dissolution into the Galactic field, an examp l e of w hich 
is the MODEST collaboration faut et all l2003t ISills et all 
2003). 

Such detailed models require stringent observational con- 
straints in the form of a precise characterization of the stellar 
content of young clusters. Investigations of the stellar popu- 
lation in young clusters hav e mostly focus ed on single stars. 
However, as pointed out by iLarsonl d200ll) . single stars only 
retain their mass from the time of formation whereas bina- 
ries retain three additional parameters, their mass ratio, angular 
momentum and eccentricity. Thus, the properties of the binary 
population can place much stronger constraints on the physical 
mechanisms underlying the star and cluster formation process. 

Ideally, one would like an accurate description of the 'pri- 
mordia l ' bina ry population. This population was defined by 
iBrownl d200ll) as "the population of binaries as established just 
after the formation process itself has finished, i.e., after the stars 
have stopped accreting gas from their surroundings, but before 
stellar or dynamical evolution have had a chance to alter the 
distribution of binary parameters". This definition is not en- 
tirely satisfactory as stellar and dynamical evolution will take 
place already during the gas-accretion phase. 

Here we revise this definition to: "the population of bina- 
ries as established just after the gas has been removed from the 
fanning system, i.e., when the stars can no longer accrete gas 
from their surroundings" . This refers to the same point in time, 
but the interpretation of the 'primordial binary population' is 
somewhat different. The term now refers to the point in time 
beyond which the freshly formed binary population is affected 
by stellar/binary evolution and stellar dynamical effects only. 
Interactions with a surrounding gaseous medium no longer take 
place. 

From the point of view of theoretical/numerical models of 
star cluster formation and evolution the primordial binary pop- 
ulation takes on the following meaning. It is the final popula- 
tion predicted by simulations of the formation of binaries and 
star clusters, and it is the initial population for simulations that 
follow the evolution of star clusters and take into account the 
details of stellar dynamics and star and binary evolution. The 
primordial binary population as defined in this pape r can be 
identifie d with the initial binary population, defined by Kroupa 
(1995b) as the binary population at the instant in time when 
the pre-main-sequence eigenevolution has ceased, and when 
dynamical evolut ion of the stellar cluster becomes effective. 
lKroupaldl995alch infers the initial binary population by the so- 
called inverse dynamical population synthesis technique. This 
method involves the evolution of simulated stellar clusters for- 
ward in time for different initial binary populations, where the 
simulations are repeated until a satisfactory fit with the present 
day binary population is found. 



Our aim is to obtain a detailed observational characteriza- 
tion of the primordial binary population as a function of stel- 
lar mass, binary parameters, and (star forming) environment. 
The most likely sites where this population can be found are 
very young (i.e., freshly exposed), low density stellar group- 
ings containing a wide spectrum of stellar masses. The youth 
of such a stellar grouping implies that stellar evolution will 
have affected the binary parameters of only a handful of the 
most massive systems. The low stellar densities guarantee that 
little dynamical evolution has taken place after the gas has 
been removed from the forming system. These constraints nat- 
urally lead to the study of the local ensemble of OB associ- 
ations. Star clusters are older and have a higher density than 
OB associations and are therefore less favorable. For example, 
in the Hyades and Pleiades, the binary population has signifi- 
cantlv changed due to dynamical and stellar evolution dKroupal 
Il995d) . Note that OB associations may s t art out as dense 
cluste rs dKroupa. Aarseth. & Hurlevl 1200 it iKroupa & Boilvl 
2002). However, they rapidly expel their gas and evolve to low- 
density systems. This halts any further dynamical evolution of 
the b inary popul a tion. 

iBrown et all dl999l) define OB associations as "young 
(550 Myr) stellar groupings of low density (< O.lM pc~ 3 ) 
— such that they are likely to be unbound — containing a 
significant population of B stars." Their projected dimensions 
range from ~ 10 to ~ 100 pc and their mass spectra cover the 
mass range from O stars all the way down to brown d warfs . 
For reviews on OB ass ociations we refer to Blaauwl dl 99 ll) 
and IBrown et all dl999l) . Thanks to the Hipparcos Catalogue 
(ESA 1997) the stellar content of the nearby OB associations 
has been established with unprecedented accuracy to a com- 
pleteness limit of V ~ 10 .5 mag, or about 1 for the stars 
in the nearest associations |de Zeeuwetal l ll999tlHoogerwerd 
2000). Beyond this limit the population of low-mass pre- 
main-sequence stars ha s been intensively studied in, e.g., the 
Sco OB2 association (Prei bisch et all Eo02t IMamaiek etall 
2002). 

The latter is also the closest and best studied of the OB 
associations in the solar vicinity and has been the most exten- 
sively surveyed for binaries. The association consists of three 
subgroups: Upper Scorpius (US, near the Ophiuchus star form- 
ing region, at a distance of 145 pc), Upper Centaurus Lupus 
(UCL, 140 pc) and Lower Centaurus Crux (LCC, 118pc). 
The ages of the subgroups range from 5 to ~ 20 Myr and 
their stellar content has been established from OB stars down 



to brown dwarfs, (for details see Ide Geu s. de Zeeuw. & Lub 
| 1989| I de Zeeuw et al l ll 999t Ide Bmiinel ll 999t iHoogerwerf 
2000J IMamaiek et ail l2002t IPreibisch et all l2002h . Surveys 
targeting the binary populatio n of Sc o OB2 include the radial- 
velocity stu dy bvlLevato et aD dl 9871) . the speckle interferome- 
try study by Kohler et all (12000 ). and the adaptive optics study 
by IShatskv & Tokovininl d2002l) . Because of their brightness 
many of the B-star members of Sco OB2 have been included 
in numerous binary star surveys, in which a variety of tech- 
niques have been employed. The literat ure data on th e binary 
population in Sco OB2 is discussed by IBrownl d200ll) and re- 
veals that between 30 and 40 per cent of the Hipparcos mem- 
bers of Sco OB2 are known to be binary or multiple systems. 
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However, these data are incomplete and suffer from severe se- 
lection effects, which, if not properly understood, will prevent 
a meaningful interpretation of the multiplicity data for this as- 
sociation in terms of the primordial binary population. The first 
problem can be addressed by additional multiplicity surveys of 
Sco OB2. In this paper we report on our adaptive optics survey 
of Sco OB2 which was aimed at surveying all the Hipparcos 
members of spectral type A and late B, using the ADONIS in- 
strument on the 3.6m telescope at ESO, La Silla. 

We begin by describing in Sect.[2]our observations, the data 
reduction procedures and how stellar components other than 
the target stars were detected in our images. These components 
have to be separated into background stars and candidate physi- 
cal companions. We describe how this was done in Sect. [3] The 
properties of the physical companions are described in Sect. [4] 
In Sect.|5]we discuss which of the physical companions are new 
by comparing our observations to data in the literature and we 
provide updated statistics of the binary population in Sco OB2. 
We summarize this work in Sect. [6] and outline the next steps 
of this study which are aimed at addressing in detail the prob- 
lem of selection biases associated with multiplicity surveys and 
subsequently characterizing the primordial binary population. 

2. Observations and data reduction 

2.1. Definition of the sample 

A census of the stellar content of the three subgroups of 
Sco OB2 based on position s, proper mo tions, and parallaxes 
was presented by Ide Zeeuw et al.l dl999h . Our sample is ex- 
tracted from their list of Hipparcos member stars and consists 
of A and B stars. In order to avoid saturating the detector during 
the observations, we were restricted to observing stars fainter 
than V ~ 6 mag. The sample therefore consists mainly of late- 
B and A stars which at the mean distance of Sco OB2 (130 pc) 
corresponds to stars with 6 mag < V < 9 mag, which translates 
to very similar limits in K$ . We observed 199 stars from the re- 
sulting Hipparcos member sample (Tabled- Not all late-B and 
A member stars could be observed due to time limitations. The 
distribution of target stars over spectral type is: 83 (157) B, 
113 (157) A, 2 (138) F, and 1 (48) G, where the numbers in 
brackets denote the total number of Sco OB2 Hipparcos mem- 
bers of the corresponding spectral type. 

2.2. Observations 

The observations were performed wi th the 
ADONIS/SHARPII+ system jBeuzit et all Il997h on the 
ESO 3.6 meter telescope at La Silla, Chile. This is an adaptive 
optics system coupled to an infrared camera with a NICMOS3 
detector array. The field of view of the 256 x 256 pixel detector 
is 12.76" x 12.76". The plate scale is 0.0495" + 0.0003" 
per pixel, and the orientation of the field is -0°.20 ± 0°.29 
(measured from North to East) 1 . Wavefront sensing was 
performed directly on the target stars. 

1 The astrometric calibrations are based on observations of the as- 
trometric reference field (6 Ori) for ADONIS/SHARPII+/NICMOS3 



10000 



1000 



I I I I I I I I 1 1 1 1 1 1 1 1 1 


■ | i i i 


: t\ 


- 


- \*\ 


- 


w 


- 


z V. X 


- 


V N 


; 


\\ \ 




V\ H 




V •. \ / 








\ •. \ 




\ -.\ 




\ *•*» 






" . V. 
w *.\ 

'* s 



0.0 0.5 1.0 1.5 2.0 2.5 
Radius (arcsec) 



Fig. 1. The radial profile of the PSF for the target 
star HIP53701. The observations are obtained using the 
ADONIS/SHARPII+ system on June 7, 2001 in J, H, and K s . 
The corresponding Strehl ratios for these observations are 18% 
in /, 26% in H and 3 1 % in K s . 

Our observing campaign was carried out in two periods: a 
period from May 3 1 to June 4, 2000, and one from June 5 to 9, 
2001. Out of the eight observing nights about 1.5 were lost due 
to a combination of bad weather and technical problems. Each 
target star is observed in the Ks (2.154 pm) band. A subset is 
additionally observed in the J (1.253 pm) and H (1.643 pm) 
band. In our search for companion stars, near-infrared observa- 
tions offer an advantage over observations in the optical. In the 
near-infrared the luminosity contrast between the primary star 
and its (often later type) companion(s) is lower, which facili- 
tates the detection of faint companions. The performance of the 
AO system is measured by the Strehl ratio, which is the ratio 
between the maximum of the ideal point spread function (PSF) 
of the system and the measured maximum of the PSF. Typical 
values for the Strehl ratio in our observations were 5-15% in 
J, 15-20% in H and 20-35% in K s (Figure^. 

Each star was observed at four complementary pointings in 
order to enhance the sensitivity of the search for close com- 
panions and to maximize the available field of view (Figure|2j. 
The two components of the known and relatively wide binaries 
HIP77315/HIP77315 and HIP80324 were observed individu- 
ally and combined afterwards. For all targets each observation 
consists of 4 sets of 30 frames (i.e. 4 data cubes). The inte- 
gration time for each frame was 200 - 2000 ms, depending 
on the brightness of the source. After each observation, thirty 
sky frames with integration times equal to those of the corre- 

(see http://www.imcce.fr/priam/adonis). Mean values for the 
period 9/1999 to 12/2001 are used. Plate scale and position angle dif- 
ferences between the two observing runs fall within the error bars. 



M.B.N. Kouwenhoven et al.: The Primordial Binary Population 



4 




Fig. 2. An example illustrating the observing strategy. 
Observations are obtained with ADONIS/SHARPII+ system 
in J, H, and K$. Four observations of each star are made, 
each observed with a different offset, in order to maximize 
the field of view. The frame size of the individual quadrants 
is 12.76" X 12.76". Combination of the quadrants results in an 
effective field of view of 19.1" x 19.1". This figure shows K$ 
band observations of the AO V star HIP75957. Two stellar com- 
ponents at angular separations of 5.6" and 9.2" (indicated with 
the black circles) are detected. These components are visible to 
the left (East) and bottom-right (South- West) of HIP75957 and 
are probably background stars. All other features are artifacts. 



sponding target star were taken in order to measure background 
emission. The sky frames were taken 10 arcmin away from the 
target star and care was taken that no bright star was present 
in the sky frame. Nevertheless, several sky frames contain faint 
background stars. The error on the target star flux determina- 
tion due to these background stars is less than 0.06% of the 
target star flux and therefore negligible for our purposes. Dark 
frames and flatfield exposures (using the internal lamp) were 
taken and standard stars were regularly observed. 

We did not use the coronograph in our setup. The corono- 
graph obscures most of the target star's flux, which compli- 
cates flux calibration of the target stars and their companions. 
Furthermore, the use of the coronograph prevents detection of 
close companions with angular separations less than 1". The 
advantage of the coronograph is that one can detect fainter ob- 
jects, which would otherwise remain undetected due to satu- 
ration of the target stars. However, the large majority of these 
faint components are likely background stars (see ^3}. 



2.3. Data reduction procedures 

The prim ary data re duction was performed with the ECLIPSE 
package (Devillard, 1997). Information about the sensitivity of 
each pixel is derived from the (dark current subtracted) flatfield 
images with different exposure times. The number of counts as 
a function of integration time is linearly fitted for each pixel. 
The results are linear gain maps (pixel sensitivity maps) for 
each filter and observing night. Bad pixel maps for each ob- 
serving run and filter are derived from the linear gain maps: 
each pixel with a deviation larger than 3cr from the median is 
flagged as a bad pixel. The average sky map is subtracted from 
the corresponding standard star and target star data cubes. The 
resulting data cubes are then corrected using the dark current 
maps, the linear gain maps, and bad pixel maps. 

Image selection and image combination was done with 
software that was written specifically for this purpose. Before 
combination, frames with low Strehl ratios with respect to the 
median were removed. Saturated frames (i.e. where the stellar 
flux exceeds the SHARPII+ linearity limit) and frames with se- 
vere wavefront correction errors were also removed. The other 
frames were combined by taking the median. Typically, about 
2 to 3 out of 30 frames were rejected before combination. 

2.4. Component detection 

We used the STARFINDER package JDiolaiti et all l2000h to 
determine the position and instrumental flux of all objects in 
the images. For each component other than the target star we 
additionally measured the correlation between the PSF of the 
target star and that of the component. Components with peak 
fluxes less than 2 to 3 times the noise in the data, or with corre- 
lation coefficients less than ~ 0.7, were considered as spurious 
detections. 

Since the PSF of the target star is generally not smooth, it 
was sometimes difficult to decide whether a speckle in the PSF 
halo of the star is a stellar component, or merely a part of the 
PSF structure. In this case four diagnostics were used to dis- 
criminate between stellar components and PSF artifacts. First, 
a comparison between the images with the target star located in 
four different quadrants was made. Objects that do not appear 
in all quadrants where they could be detected were considered 
artifacts. Second, a similar comparison was done with the indi- 
vidual uncombined raw data frames. Third, a radial profile was 
fit to the PSF and subtracted from the image to increase the 
contrast, so that the stellar component becomes more obvious. 
And finally, a comparison between the PSF of the target star 
and the PSF of the star previous or next in the program with a 
similar position on the sky was made by blinking the two im- 
ages. These two PSFs are expected to be similar and therefore 
faint close companion stars can be detected. 

In total we detect 151 stellar components other than the tar- 
get stars. A significant fraction of these 151 components are 
background stars (see §[3}- All other components are compan- 
ion stars. Figure [3] shows the K$ magnitude of the compan- 
ion and background stars as a function of angular separation 
p between the target star and the companion or background 

Stars. The Same plot also shows AK S = ^.component - ^.primary 
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as a function of p. Previously unknown companion stars (see 
§ 15 - 1 1 and known companions are indicated with the dots and 
the plusses, respectively. The background stars (see § are 
represented with triangles. This figure shows that the detec- 
tion probability of companion stars increases with increasing 
angular separation and decreasing magnitude difference, as ex- 
pected. The lower magnitude limit for component detection is 
determined by the background noise; the lower limit for the an- 
gular separation is set by the PSF delivered by the AO system. 
Note that not all detected components are necessarily compan- 
ion stars. A simple criterion is used to separate the companion 
stars (Ks < 12 mag) and the background stars (Ks > 12 mag). 
See §|3]for the motivation of this choice. 

The estimated detection limit as a function of angular sep- 
aration and companion star flux is represented in Figure[3]with 
a solid curve. The detection limit is based on Ks band mea- 
surements of simulated companion stars around the B9 V tar- 
get star HIP65178. In the image we artificially added a second 
component, by scaling and shifting a copy of the original PSF, 
thus varying the desired flux and angular separation of the sim- 
ulated companion. We sampled 30 values for the companion 
flux linearly in the range 9.35 mag < Ks < 15.11 mag, and 40 
values for the angular separation between primary and compan- 
ion linearly in the range 0.05" < p < 2". For each value of the 
angular separation we determined the minimum flux for which 
the companion star is detectable (i.e. with a peak flux of 2 - 3 
times the noise in the data). In order to minimize biases, we re- 
peated this procedure for four different position angles and av- 
eraged the results. The standard deviation of the detection limit 
is 0.07" at a given magnitude and 0.35 mag at a given angu- 
lar separation. Similar simulations with other target stars (e.g., 
HIP58859, Strehl ratio = 8%; HIP76048, Strehl ratio = 27%) 
show that the detection limit for the observations with different 
Strehl ratios is comparable to that of HIP65178 (Strehl ratio = 
40%). All target stars are of spectral type A or late-B and have 
comparable distances. The detection limit shown in Figure[3]is 
therefore representative for our sample. 

2.5. Photometry 

Near-infrared magnitu des for the observed standa r d star s 
are taken from [y an der Bliek. Manfroid. & Bouchet ( 1996), 
ICarter & Meadows! dl995h . and the 2MASS catalog. All 
magnitudes are converted to the 2MAS S syste m using the 
transformation formulae of Carpenter (2001). The stan- 
dard stars HD101452, HD190285 and HD96654 were found 
to be double. The fluxes for the primary and secondary 
of the standard stars HD101452 and HD96654 are added 
for calibratio n, since the second a ries a re unresolved in 
Ivan der Bliek. Manfroid. & Bouchetl Jl996l) . For HD190285, 
only the flux of the primary is used for calibration since the 
secondary is resolved in the 2MASS catalog. 

The zeropoint magnitudes are calculated for each observing 
night and each filter individually. There are non-negligible pho- 
tometric zeropoint differences between the four detector quad- 
rants. Application of the reduction procedure to simulated ob- 
servations (but using the ADONIS/SHARPII+ dark frames and 



Hatfield exposures) also shows this zeropoint offset, unless the 
flatfielding was omitted. Hence, the zeropoint offset between 
the four quadrants is caused by large-scale variations in the flat- 
field exposures. The four detector quadrants were consequently 
calibrated independently. 

Since the distribution of observed standard stars over air- 
mass is not sufficient to solve for the extinction coefficient k, 
we use the mean extinction coefficients for La Silla 2 : kj = 
0.081 mag/airmass, kn = 0.058 mag/airmass and kfc s = 
0.113 mag/airmass. The standard stars have similar spectral 
type as the target stars, which allows us to neglect the color 
term in the photometric solution. 

3. Background stars 

Not all of the 151 stellar components that we find around the 
target stars are companion stars. The chance that a foreground 
or background star is present in the field is not negligible. There 
are several techniques that can be applied to determine whether 
a detected component is a background star or a probable com- 
panion star. 

Separation of the objects based on the position of the com- 
ponent in the color-magnitude diagram (CMD) is often ap- 
plied to discriminate backgr ound stars from companion stars 
(see, e.g., IShatskv & Tokovinirll2002l) . Given a certain age of 
the population, a companion star should be found close to the 
isochrone that corresponds to that age. When age and distance 
spread are properly taken into account, this method provides an 
accurate estimate of the status of the object. We are unable to 
apply this technique to our data since only a small fraction of 
our observations is multi-color (see § !4.2l and FigurefTol. 

Following IShatskv & Tokovininl J2002I) we use a simple 
brightness criterion to separate the companion stars (Kg < 
12 mag) and the background stars {Ks > 12 mag). At the 
typical distance of a Sco OB2, an M5 main sequence star 
has approximately Ks = 12 mag. The three subgroups of 
Sco OB2 are located at distances of 145, 140, and 1 18 pc. We 
assume an age of 5 Myr for US jde Geus. de Zeeuw. & Lubl 
ll989tlPreibisch & Zinneckelll999l) . and 20 Mvr for UCL and 
LCC llMamaiek et all 120021) . Using the isochrones described 
in § 14.31 we find that stars with Ks = 12 mag correspond to 
stars with masses of 0.08 M o ,0.18 M Q , and 0.13 M for US, 
UCL, and LCC, respectively. Using the Ks = 12 mag criterion 
we conclude that 77 out of 151 components are likely back- 
ground stars. The remaining 74 components are thus candidate 
companion stars. 

Although this Kg = 12 mag criterion works well for 
stars, it is not very accurate below the hydrogen burning limit. 
Mar tin et alj j2004) found 28 candidate brown dwarf members 
of the US subgroup using /, J, and K photometry from the 
DENIS database. All members have 10.9 mag < K < 13.7 mag 
and their spectral types are in the range M5.5-M9. According 
to the Ks = 12 mag criterion, 19 of these brown dwarfs would 
be classified as companion stars and 8 as background stars (one 
brown dwarf does not have a K magnitude entry). A significant 

2 see http : //www . Is . eso . org/lasilla/sciops/ntt/sof i 
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Fig. 7. The distribution of 74 detected candidate companions (top) and the 77 presumed background stars (bottom) over Galactic 
coordinates. The number of companions and background stars found per target star is indicated by the size of the diamonds and 
triangles, respectively. Hipparcos member stars that were not observed are shown as dots. The borders of the three subgroups 
of Sco OB2 (Upper Scorpius, Upper Centaurus Lupus, and Lower Centaurus Crux) are indicated by the dashed lines. The 
background stars are concentrated towards the location of the Galactic plane (see also Figure[8}. 



fraction of possible brown dwarf companion stars to our ob- 
served target stars would be classified as background stars. 

In this study we cannot determine with absolute certainty 
whether a detected component is a companion star or a back- 
ground star. For example, a foreground star with K$ > 12 mag 
will not be excluded by the criterion described above. However, 
the differences between the companion star distribution and 
the background star distribution can be used to find out if the 



results are plausible or not. For this purpose, we derive the 
expected distribution of background stars per unit of angular 
separation P(p), position angle Q(tp), and Ks band magnitude 
R(K$) for the background stars. 



Background stars are expected to be uniformly distributed 
over the image field with a surface density of I stars per unit 
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of solid angle. The number of observed background stars P(p) 
per unit of angular separation p is given by 



(1) 



where Q p is the part of the field of view contained within a 
circle with radius p centered on the reference object (the tar- 
get star). In calculating the expected P(p) for our observations, 
we assume that the reference object is located in the center of 
the field of view. For a square field of view with dimension L, 
Equation^then becomes 



P(p) oc . 



2np for < p 

8p(| - arccos j-) f° r < P 
for p > 



< L/2 

< L/V2 
L/V2 



(2) 



For our observing strategy (FigurefJJl we haveL = §■ 12.76" = 
19.14". For p < L/2 = 9.6" the background star density for 
constant K$ is proportional to p, after which it decreases down 
to at p — L/y/2 = 13.6". Here we have made the assumption 
that each target is exactly in the center of each quadrant. 

We can also calculate Q(tp), the expected distribution of 
background stars over position angle tp. For p < L/2, Q(tp) 
is expected to be random. For L/2 < p < L/V2 one expects 
Q(<p) = L 2 /8 cos 2 tp for 0° < <p < 45° 3 . Q(ip) is undefined 
for p > L/y/2. The above calculations apply to an ideal situa- 
tion of our observing strategy. In practice, the target star is not 
always exactly centered in one detector quadrant, which influ- 
en ces P(p) and Q(tp) for p ^ Lh j2. 

LShatskv & Tokovininl ll2002) analysed the background stars 
in their sky frames taken next to their science targets in 
Sco OB2. We use their data to derive R(Ks), the expected 
K s magnitude distribution for background stars. We fit a sec- 
ond order polynomial to the logarithm of the cumulative K s 
as a function of Ks of th e background stars (Figure 4 in 
IShatskv & Tokovininl Eoo3> . In our fit we only include back- 
ground stars with 12 mag < Ks < 16 mag since that is the Ks 
magnitude range of the presumed background stars. The Ks 
magnitude distribution f(Ks ) is then given by the derivative of 
the cumulative Ks distribution. 

Figure |4] shows the companion star and background star 
distribution as a function of angular separation. The two dis- 
tributions are clearly different. The companion stars are more 
centrally concentrated than the background stars. The solid 
curve represents Pip), normalized in such a way that the num- 
ber of stars corresponds to the expected number of back- 
ground stars with 12 mag < Ks < 16 mag (see below and 
Figure |5}. More background stars are present at angular sepa- 
ration 4" < p < 6" relative to what is expected. This is con- 
firmed with a Kolmogorov-Smirnov (KS) test. The KS test is 
usually quantified with the "KS significance", which ranges 
from to 1 and corresponds to the probability that the two 
datasets are drawn from the same underlying distribution. For 
this KS comparison we only consider those background stars 



3 This formula is also valid for 90° < tp < 135°, 180° < tp < 125°, 
and 270° < ip < 315°. Q(tp) = L 2 /8 sin 2 tp for all other position angles 
between 0° and 360°. 



with p < L/V2 = 9.6". For background stars with p > 9.6" 
we expect the predicted distribution to differ slightly from the 
observed distribution (see above). We find a KS significance 
level of 1 .8 x 1CT 4 , which implies that the excess of background 
stars with 4" < p < 6" is real, assuming that the simple model 
for the angular separation distribution is correct. The excess of 
close background stars might be caused by faint (Ks > 12 mag) 
companion stars which are misclassified as background stars. 

The position angle distribution of companion stars and 
background stars is shown in Figure |6] The objects are ex- 
pected to be randomly distributed over position angle for p < 
L/2 = 9.6", which is illustrated in Figure |6] For objects with 
p < L/2 = 9.6" we find that, as expected, the position an- 
gle distribution is random, for both the companion stars (KS 
significance level 0.86) and the background stars (KS signif- 
icance level 0.84). For angular separations larger than 9.6" 
the objects are primarily found in the 'corners' of the image 
(tp = 45°, 135°, 225°, 315°). 

The expected number of background stars per unit of an- 
gular separation and unit of magnitude, R(Ks)P(p) is plotted 
in Figure [5] We observe 19 background stars in the region 
14 mag < Ks < 15 mag and 5" < p < 13" and normalize 
the theoretical background star density accordingly. The bright- 
ness of the PSF halo can explain the smaller number of back- 
ground stars than expected for p ^ 2.5" and Ks ^ 15 mag. 
The predicted strong decline in background stars for p ^ 10" 
is present in the observations. Few background objects with 
12 mag < Ks < 14 mag are expected in the angular separation 
range of 2" -4". 

The distribution of the detected candidate companions and 
background stars on the sky is shown in Figure The cov- 
ered fields close to the Galactic plane include a larger number 
of background stars, as expected. This is even better visible in 
Figure[8] which shows the distribution of observed objects over 
Galactic latitude. The distribution of Sco OB2 Hipparcos mem- 
ber stars over Galactic latitude is similar to that of the compan- 
ion stars. 

In principle, it is also possible that another (undocumented) 
member star of Sco OB2 projects close to the target star. The 
status of these stars (companion star or "background star") can- 
not be determined w ith the CMD or the Ks = 12 mag criterion. 
IPreibisch et alJ J2002I) estimate the number of member stars in 
the US subgroup in the range 0.1 - 20 M Q to be about 2500. 
This corresponds roughly to a stellar surface density of 8 x 10~ 7 
stararcsec~ 2 . The chance of finding another US member star 
(i.e. not a binary companion) in a randomly pointed observa- 
tion of US is of the order 0.03%. The detection probabilities 
for the other two subgroups of Sco OB2 are similarly small 
and negligible for our purpose. 

4. Properties of the companion stars 

4.1. General properties 

The angular separation distribution of the companion stars is 
centrally concentrated (Figures |4] and |6j. The error in the an- 
gular separation is typically 0.0015". Faint (Ks ^ 11 mag) 
objects with angular separation less than p * 0.8" cannot be 
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detected in the presence of the PSF halo of the bright pri- 
mary star, while the upper limit of 13.5" is determined by the 
field of view. The distribution over position angle is random 
for p < L/2 = 9.6", as expected (see §[3}. The median error 
in position angle is 0.03°, increasing with decreasing angular 
separation. 

Two observed systems are clearly hierarchical triples 
where a double 'secondary' orbits the primary star. These are 
HIP68532, where the double 'secondary' has p = 3.1" and 
ip = 290°, and HIP69113, with <p = 65° and p = 5.4". 
HIP69113 has a third companion star with p = 28.6" and 
tf = 35° iDuflot. Fieon. & Mevssonnierl dl995h . This star is 
not observed in our survey because of the large distance from 
the primary star. The nature of the H IP69113 system is ex- 
tensively discussed by Huelamo et al. (2001). Figure [9] shows 
HIP68532 and HIP69113 with companion stars. The target 
stars HIP52357, HIP61796, and HIP81972 also have two com- 
panion stars. The two companion stars of each of these systems 
have approximately the same angular separation and position 
angle, but it is unclear whether these systems are hierarchical 
in the way described above. 

The 199 observed primaries have a Ks magnitude range 
of 5.2 - 8.3 mag. The Ks magnitudes of the companion stars 
(Figure range from Ks = 6.4 mag to 12.0 mag, the upper 
limit resulting from the criterion that was used to remove the 
background stars from the sample. The primaries in our sample 
all have similar spectral type, distance, and interstellar extinc- 
tion. The magnitude of the companions relative to that of their 
primaries, AKs, spans the range from 0.0 mag to 6.0 mag. The 
median error is 0.05 mag in J, 0.04 mag in H, 0.07 mag in Ks, 
and 0.075 mag in / - K. 

Properties of the 199 target stars, 74 candidate companion 
stars, and 77 presumed background stars are listed in Tabled 
The name of the star is followed by the J, H, and Kg mag- 
nitude. If the magnitude is derived from measurements done 
under non-photometric conditions, a remark is placed in the 
last column. The spectral type of each primary is taken from 
the Hipparcos catalogue. The angular separation and position 
angle (measured from North to East) are derived from the com- 
bination of all available observations for a particular star. For 
each star we also list the status (p = primary star, c = candi- 
date companion star, nc = new candidate companion star, b = 
background star) and the subgroup membership. 

4.2. Color-magnitude diagram and isochrones 

In § [5]we mentioned that it is possible to separate background 
stars and companion stars using their position in the CMD. In 
this section we derive the absolute magnitude Mk s for all ob- 
served objects. We additionally derive colors for objects with 
multi-color observations. Only measurements obtained under 
photometric conditions are used. We construct a CMD and de- 
termine whether the observed components are background stars 
or companion stars by comparing their position in the CMD to 
that of the isochrones. In § 14. 31 we will use the absolute mag- 
nitude Mk s and the isochrones to derive the masses of the pri- 
mary and companion stars. 



We calculate the absolute magnitude Mk s using the dis- 
ta nce D and extinction Ak s for each target star individually. 
In lde Bruiinel (fl999) the distance to the Sco OB2 member stars 
is derived from the secular parallax 7r sec . The secular parallax 
is calculated from the observed positions and proper motions 
of member stars, where the fact that all stars in an association 
share the same space motion is exploited. The secular paral- 
lax can be up to more than two times as precise compared to 
the Hipparcos parallax. Secular parallaxes with g — 9 are used 
when available, and with g = oo o therwise (g mea sures the 
model-observation discrepancy; see Ide Bruiind. fl999l for de- 
tails}. The interstellar extinction is also taken from Ide Bruiind 
dl999h. Fo r the s tars that do not have an extinction entry in 
Ide Bruiind (Ep99), an estimate of the extinction is calculated 
using: 

A Ks =R v x E(B - V)/9.3, (3) 

whereby = 3 .2 is the standard ratio of total-to-selective extinc- 
tion (see, e.g., ISavage & Mathislll979l) . The value 9.3 is the ra- 
tio between V band and Ks band extinction jMathisl fl990). and 
E(B - V) = ( B-V)-(B-V)qis the color excess. (B - V) is the 
color listed in de Bruiinel dl999l) . The theoretical color (B-V)o 
is derived from the spe ctral type using the bro adband data for 
main sequence stars in Kenyon & Hartmann (1995). The ex- 
tinction for the UCL member HIP68958 (spectral type Ap...) is 
not listed in de Bruiine ( 1999) and could not be derived from 
Kenvon & Hartmann ( 1995). We therefore take Avhip68958 = 
0.063 mag. This is the median Ay value fo r those UCL mem - 
ber stars which have full photometric data in de Bruiine (1999). 

Only for a subset of observed targets observations in more 
than one filter are available; the corresponding CMDs are plot- 
ted in Figure [K)] The observed primary stars are all of simi- 
lar spectral type but show scatter in Mk s , J — K, and H - K. 
This can be explained by the errors in near-infrared photom- 
etry, which are typically 0.07 mag for our observations. The 
uncertainties in the parallax and reddening introduce an addi- 
tional error when calculating the absolute magnitudes. Median 
errors in M; and Mk are consequently about 0.36 mag. The 
colors J - K and H - K have median errors of 0.02 mag and 
0.06 mag, respectively. 

Given the age of the three Sco OB2 subgroups, not all 
objects are expected to be positioned on the main sequence. 
For the young age of Sco OB2, stars of spectral type G or 
later have not reached the main sequence yet, while the O stars 
have already evolved away from the main sequence. The age 
differences between the Sco OB2 subgroups are relatively 
well-determined. Different values for the absolute ages are 



derived from the kinematics of the subgroups ||B. 


aauw 


19641 Il978h and stellar evolution dde Zeeuw & BrandT 


1985 


Ide Geus. de Zeeuw. & Lubl Il989t Ipreibisch et all 


2002 


Mamaiek et al].|2Q02: Sartori. Lepine. & Diasl 12003). The ase 



of the US subgroup i s 5 Myr, without a significant age spread 
dde Zeeuw & Brandt ll985Ude Geus. de Zeeuw. & Lubl Il989t 
Preibisch et all 120021) . Mamaie k et aid d2002) derive an age 
between 15 and 22 Myr for members of UCL and 17 and 
23 Myr for members of LCC. In this paper we adopt an age of 
20 Myr for both UCL and LCC. 
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Table 1. Results from the ADONIS adaptive optics survey among 199 Sco OB2 member stars. The columns list: Hipparcos 
number, observed J, H and Ks magnitude, spectral type of the primary, angular separation (arcsec), position angle (degrees; 
measured from North to East), and status (p = primary; c = companion star; nc = new companion star; b = background star). 
The last two columns list the subgroup of which the system is a member (US = Upper Scorpius; UCL = Upper Centaurus Lupus; 
LCC = Lower Centaurus Crux), and a flag in case the observing conditions were not photometric. The Hipparcos member 
stars HIP77315 and HIP77317 are a common proper motion pair (see the Double Star Catalog). In this paper the two stars are 
treated as separate primaries and not as one system. Table 1 is presented in its entirety in the electronic edition of Astronomy & 
Astrophysics. A portion is shown here for guidance regarding its form and content. 



We construct isochron es whi ch are very similar to those de- 
scribed in Preibisch et alJ 1120021) . For the low-m ass stars (M < 
0.7 Mp ) we use exactly the sam e isochrone as [Preibisch et alJ 
(2002J): we use the models from LC habrier et al. (2000), which 
are based on those from Ba raffe et alJ ( 119981) For 0.02 < 
M/M < 0.7 we take the models with [M/H] = and 
£mix = Hp, where L mlx is the mixing length and Hp the pres- 
sure scale height. For 0.7 < M/M < 1 we use the models with 
[ M/Hl = and Lrni* = \.9Hp. We use the models described 
in lPalla & StahleJ dl999h for 1 < M/M < 2. Near-infrared 
magnitude tracks, derived using the procedure described in 
iTesti. Palla. & Nattal dl998h . were kindly provide d by F. Palla 
Finally , for M/M > 2 we use t he mo dels from lGirardi et alJ 
J2002I) (based on iBertelli et"ail 1 1994 with Y = 0.352 and 
Z = 0.05, where Y and Z are the helium and metal abundance, 
respectively. 

The resulting 5 Myr isochrone (for US) and 20 Myr 
isochrone (for UCL and LCC) are plotted in Figure ^3 Both 
isochrones are continuous. Figure^2 snows the relationship be- 
tween stellar mass and ^-band absolute magnitude for 5 Myr 
and 20 Myr populations. The gray-shaded area in Figures fTol 
andll His enclosed by the 15 Myr and 23 Myr isochrones, and 
shows the effect of the age uncertainty in the UCL and LCC 
subgroups. 

The positions of the primary and companion stars in the 
CMD are consistent with the isochrones (except the companion 
of the LCC member star HIP53701). No obvious correlation of 
the position of these primaries in the CMD with spectral type 
is present due to the errors in the colors and magnitudes. The 
presumed background stars are located far from the isochrones 
in the CMD, supporting our hypothesis that they are indeed 
background stars. 



In Figure ^3 we added the can didate brown dwarf mem- 
bers in US from Martin et al. (2004). We estimate the absolute 
magnitude of these brown dwarfs using the mean US distance 
of 145 pc and the mean extinction Ak = 0.050 mag (see be- 
low). The colors and magnitudes of these objects are consistent 
with the 5 Myr isochrone (solid curve), as is expected for US 
member stars. 

4.3. Masses and mass ratios 

We derive the mass of the primary and companion stars from 
their Mk s magnitude using the isochrones described in § 14.21 
( Fi gures fTol and [Til . We use the 5 Myr isochrone for US and 
the 20 Myr isochrone for UCL and LCC. The conversion from 
Mk s to mass is most accurate for B- and A-type stars, where 
the absolute value of the derivative dM/dM^ is small. We 
find masses between 1.4 M and 7.7 M Q for the primaries and 
masses between 0.1 M and 3.0 M for the companions. Given 
the error AMk s on Mk s , we calculate for each star the masses 
M_ and M+ that correspond to M Ks - AM Ks and M Ks + AM Ks . 
We define the error on the mass to be j(M+ - MJ). The median 
of the error in the mass, which is a lower limit for the real error, 
is 0.4 M for the primaries and 0.1 M for the companions. 

Figure^l(top) shows the mass distribution of the compan- 
ion stars and of the primary stars to which they belong. The 
mass distribution of the observed target stars without compan- 
ions is similar to the latter. The plot showing companion star 
mass as a function of angular separation (Figure bottom) 
closely resembles Figure[3] The reason for this is that the com- 
panion star mass is derived from the Ks magnitude. 

The mass ratio q and error are calculated for each primary- 
companion pair. The median formal error on the mass ratio 
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Fig. 10. The color-magnitude diagram of the primaries, companion stars and background stars for which we obtained 
ADONIS/SHARPII+ multi-color AO observations. The Mk s magnitudes are corrected for distance and extinction. Companion 
stars in the three different subgroups are indicated with plusses (US), asterisks (UCL) and crosses (LCC). Target stars and back- 
ground stars are indicated with diamonds and triangles, respectively (see §[3j. The curves represent isochrones of 5 Myr (solid 
curve) and 20 Myr (dotted curve). The 15 Myr and 23 Myr isochrones enclose the gray-shaded area and represent the uncertainty 
in the age of the UCL and LCC subgroups. The mass scale is in dicated for the 20 Myr isochrone (left) and the 5 Myr isochrone 
(right). The brown dwarf candidates identified bv lMartm et alJ d2004l) are indicated as dots. The median formal errors are indi- 
cated as error bars in the top-right corner of each plot. The locations of the target stars and companion stars are consistent with 
the isochrones (within the error bars). The presumed background stars are located far away from the isochrones, implying that 
our criterion to separate companion stars and background stars is accurate. 



is 0.07 and is to first order proportional to that of the com- 
panion star mass. The real error on the mass ratio is most 
lik ely larger than the formal e rror stated above. As indicated 
by Shatskv & Toko vininl J2002I) . the mass ratio of a system pri- 
marily depends on the magnitude difference between primary 
and secondary. 

The mass ratio distribution is shown in Fi gures fL3*l and ITU 
The smallest mass ratios are found for the companions of 
HIP80474 (q = 0.022 ± 0.008) and HIP77911 (q = 0.035 ± 
0.012). This corresponds well with the expected minimum 
value for q: a typical primary in our sample has spectral type 
AO V (2 M ), while the least massive companions considered 
have Ks = 12 mag (~ 0.1 M ), giving q = 0.05. The systems 
with large mass ratio (q > 0.9) are HIP80324 (q = 0.91+0.23), 
HIP50520 (q = 0.93 ± 0.19), HIP80238 (q = 0.94 ± 0.29), 
HIP64515 (q = 0.94 ± 0.15), HIP61639 (q = 0.95 ± 0.15), and 
HIP52357 (q = 1.00 ± 0.34). Note that for these systems the 
possibility that q > 1 is included. Due to the error in the mass 
it is not possible to say which star in these systems is more 
massive, and hence, which is the primary star. Although the 
uncertainty in the age of UCL and LCC is not negligible, the 



effect on the mass ratio distribution is small (gray-shaded areas 
in Figure fT4li . 



Following Shat skv & Tokovininl (2002) we fit a power-law 
of the form f(q) = q~ r to the mass ratio distribution. We find 
that models with F = 0.33 fit our observations best (KS signif- 
icance level 0.33). The error in F due to the age uncertainties 
in the UC L and LCC subgroups i s 0.02 . This is in good agree- 
ment with lShatskv & Tokovininl ll2002l) . who surveyed 1 15 B- 
type stars in Sco OB2 for binarity and observed 37 physical 
companions around these stars. They find F = 0.3 - 0.5 for 
their mass ratio distribution. We observe an excess of systems 
with q ~ 0. 1 with respect to what is expected for models with 
T = 0.33. This could partially be explained by bright back- 
ground stars that are misclassified as companion stars, but this 
is unlikely to be an important effect (see §[3}. The observational 
biases (e.g. the detection limit) are not taken into account here. 
For a more detailed description of the effect o f observational bi- 
ase s on the mas s ratio distribution we refer to Hogeveen ( 1990) 
and lToutllll99ll) . 
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The cumulative mass ratio distribution F(q) can also be de- 
scribed as a curve consisting of two line segments. We fit and 
find that the mass ratio distribution closely follows the function 

f 2.63^ - 0.067 for 0.03 < q < 0.19 
{q) ~ \ Q.llq + 0.294 for 0.19 < q < 0.97, ' 

with a root-mean-square residual of 0.019 (see FigurelTH. 

We investigate whether the observed mass ratio distribution 
could be the result of random pairing between primary stars 
and companion stars, su ch as observed for solar-typ e stars in 
the solar neighbourhood (Duauennov & M avorl.fl 99 lh . To this 
end, we use Monte Carlo simulations to calculate the mass ratio 
distribution that is expected for random pairing. 

The current knowledge about the brown dwarf pop- 
ulati on in Sco OB2 is inc omplete (e.g., Table 2 in 
IPreibisch. Stanke. & Zinneckerl 120031) . Therefore we make a 
comparison with the mass ratio distribution resulting from two 
different initial mass functions (IMFs) (IMF_o.3 and IMF2.5), 
which differ in slope dN /dM for substellar masses. 



IMF„ 



dN 
dM 



' M- 


' 


for 


0.02 <M/M 


< 0.08 


M~ 


-0.9 


for 


0.08 < M/M Q 


< 0.6 


M~ 


■2.8 


for 


0.6 < M/M Q 


< 2 


Mr 


•2.6 


for 


2 < M/M 


< 20, 



(5) 



For M > 0.10 M n IMF„ is equa l to t he IMF in US that was 
derived by IPreibisch et all J2002). The IPreibisch et ail (§002) 
IMF is extrapolated down to M = 0.08 M . For 0.02 < 
M/ Mp < 0.08 we consider a _ -0.3 jKroupalHool) and a = 
2.5 ( IPreibisch. Stanke & Zinneckerl 120031 Fit T). Mo s t other 
models described in Preibisch. Stan ke. & Zinneckerl (|2003) 
have -0.3 < a < 2.5 for 0.02 < M/M Q < 0.08. 

We draw the primary and secondary stars from IMF„. 
Primary and secondary masses are independently drawn from 
the mass range [0.02, 20] M and are swapped if necessary, so 
that the primary star is the most massive. Only those systems 
with primary mass larger than 1 .4 M and smaller than 7.7 M 
are considered. We calculate the resulting mass ratio distri- 
bution. A KS comparison between the mass ratio distribution 
resulting from random pairing and the observed distribution 
shows that the random-pairing hypothesis can be rejected for 
both IMF_o.3 and IMF2.5. This is clearly visible in Figurell4l 

Random pairing requires more systems with low mass ratio 
than we observe. The error in the mass ratio (ranging between 
0.008 and 0.34) cannot compensate for the lack of faint com- 
panions. Even if all close (p < 4) background stars are treated 
as low-mass companions, random pairing is excluded. Hence, 
there is a deficit of low-mass companions around late-B and A 
type stars compared to what is expected from random pairing. 

Several observed primary stars are actually unresolved bi- 
nary systems. Since we derive the mass for this 'primary' using 
the total Ks flux of such an unresolved system, we introduce an 
error in the mass and hence the mass ratio. Out of the 199 ob- 
served primary stars 18 are known spectroscopic binaries (14 
in US, 3 in UCL, and 1 in LCC). These binary systems are too 
tight to be resolved in our AO survey. We observe companion 
stars around three of these primaries: HIP77911, HIP77939, 
and HIP79739 (all in US). For these systems we expect the 



mass ratio, which we define here as the mass of the companion 
star divided by the total mass of the inner part of the system, 
to deviate significanly from the value derived from our near- 
infrared observations. 

A closer inspection of Figure [3] reveals a lack of stellar 
components with p < 3.75" and 12 mag S5 Ks ^ 14 mag. 
If companion stars or background stars with these properties 
exist, we should have seen them. In §[3] we showed that few 
background stars are expected in this r egion. For example, the 
brown dwarfs (M5.5-M9) detected bv iMartm et alJ (12004 in 
US have 11 mag H K S5 14 mag. If similar objects within 
p < 3.75" were present, they would have been detected. This 
'gap' is not the result of observational biases: the detection 
limit shows that objects with these properties are detectable. 

In order to find out if this 'gap' is expected, we gener- 
ate a stellar population using a Monte Carlo process and two- 
dimensional two-sample KS tests. For comparison with the 
simulated data we only consider the 59 observed companion 
and background stars with p < A" . Objects with p > 4" are 
not relevant for this test and might introduce biases due to ob- 
servational selection effects. We randomly draw primary and 
secondary masses in the mass range [0.02, 20] M . We select 
only those binaries for which the primary mass is in the range 
[1 .4, 7.7] M . Sets of 59 binary systems are created for IMF_o.3 
and IMF2.5 (Equation |3 ■ The masses are converted into K$ 
magnitudes using a distance of 145 pc and the 5 Myr isochrone 
described in § 14.21 The angular separations for the secondaries 
are drawn from a uniform distribution in [0",4"]. We create 
10 3 realizations which are compared to the observations. For 
IMF_o.3 we find a mean KS significance level of 3.5 x 10~ 3 , 
while for IMF2.5 we find 3.2 x 10~ 3 . We repeated the procedure 
described above with another angular separation distribution: 
dN/dp ocpforpe [0",4"] (Opik's law). We find mean KS 
sigificance levels of 8.2 x 10~ 5 and 6.5 x 10~ 5 for IMF_ .3 and 
IMF2.5, respectively. It is therefore unlikely that the observed 
distribution (including the 'gap') is drawn from the IMFs and 
angular separation distributions as described above. 

Several faint close components with p < 3.75" below the 
'gap' (Ks ^ 14 mag) are detected. These objects are found next 
to the target stars HIP61265, HIP67260, HIP73937, HIP78968, 
HIP79098, HIP79410, and HIP81949. The Strehl ratios in K s 
of the corresponding observations (~ 20%) are typical. The 
detection of these close, faint components therefore cannot be 
expained by a better performance of the AO system during the 
observation of these objects. The target stars corresponding to 
these objects are not particularly faint; the luminosity contrast 
between the target star and the faint close components is typ- 
ical. These objects are background stars according to our se- 
lection criterion. However, there is a possibility that these ob- 
jects are close brown dwarfs with masses ^ 0.05 M , where the 
conversion from Ks to mass is strongly dependent on the age. 
This result implies that A and late-B stars do not have close 
companions with a mass less than about 0.1 M , unless the 
assumed background stars are physical companions (and thus 
brown dwarfs). If so, a gap would be present in the companion 
mass distribution. 
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Reference 



Detection method 



Alencaretal. (2003) 

Balega et al. (19941 

Barbi er-Brossat. Petit. & Figon (1994) 

Batten. Fletcher. & MacCarthv (1997) 

Couteau (19951 

The Double Star Library 

Duflot. Fiaon. & Mevssonnier (1995) 

HartkoDf. McAlister, & Mason (2001) 

Jordi et al. (1997) 

The Hipparcos and Tycho Catalogues 

Kraicheva et al. (19891 

Malkov (19931 

Mason (19951 

McAlister et al. (19931 

Miscellaneous, e.g. SIMBAD 

Miuraetal. (19921 

Pedoussa ut et al. (19961 

Shats kv & Tokovinin (20021 

Sowell & Wilson (19931 

Svechnikov & Bessonova (1984) 

Tokovinin (19971 

Worlev & Douglass (19971 



Spectroscopic 
Visual 

Spectroscopic 

Spectroscopic 

Combination 

Combination 

Spectroscopic 

Visual 

Eclipsing 

Astrometric 

Spectroscopic 

Combination 

Visual 

Visual 

Combination 
Visual 

Spectroscopic 

Visual 

Visual 

Combination 
Combination 
Combination 



Table 2. References to literature data with spectroscopic, astro- 
metric, eclipsing, and visual binaries in Sco OB2. Using these 
data we find that 33 out of the 74 candidate companion stars 
in our dataset were already documented in literature, while 41 
were previously unknown. 



Type 


B1-B3 


B4-B9 


A 


F 


GKM 


Upper Scorpius 


Single stars 


3 


13 


20 


17 


10 


Visual 


15 


9 


12 


3 


2 


Astrometric 





3 


1 


2 


3 


Spectroscopic 


11 


12 











Upper Centaurus Lupus 


Single stars 


3 


22 


41 


43 


22 


Visual 


15 


23 


29 


7 


8 


Astrometric 


2 


1 


3 


5 


2 


Spectroscopic 


11 


6 


1 


1 





Lower Centaurus Crux 


Single stars 


3 


16 


31 


47 


14 


Visual 


7 


13 


24 


10 


4 


Astrometric 


1 


3 


3 


4 


4 


Spectroscopic 


3 


3 


1 








Scorpius OB2 


Single stars 


9 


51 


92 


107 


46 


Visual 


37 


45 


65 


20 


14 


Astrometric 


3 


7 


7 


11 


9 


Spectroscopic 


25 


21 


2 


1 






Table 3. Binarity data for Sco OB2, showing the number of 
companion stars found with the different techniques. Literature 
data and the new companion stars described in this paper are 
included. The terms 'visual', 'astrometric', and 'spectroscopic' 
pertain only to available observations and not to in intrinsic 
properties of the systems. Note that each companion of a mul- 
tiple system is included individually. 



5. Comparison with literature data 



5.1. New companions 

We compiled a list of known companion stars to all Hipparcos 
members of Sco OB2 using literature data on binarity and mul- 
tiplicity. This includes spectroscopic, astrometric, eclipsing, 
and visual binaries. The references used are listed in Table|2] 

Since the angular separation and position angle of compan- 
ion stars change as a function of time, one has to take care that 
a 'newly' detected companion is not a displaced known com- 
panion star. It is therefore interesting to estimate how fast the 
angular separation and position angle change for the observed 
companion star. We obtained the primary mass M p and com- 
panion star mass M s in (§ I4.3> thus we can estimate the orbital 
period of the companion stars using Kepler's third law. If the 
orbit is circular and the system is seen face-on, the orbital pe- 
riod is given by P = tJ4tt 2 {Dp) 3 /G{M p + M s ), where D is the 
distance to the system, p the angular separation between pri- 
mary and companion star, and G is the gravitational constant. 
In general, orbits are eccentric and inclined. However, the effect 
of nonzero eccentricity and inclination on the period is most 
likely only of the order of a few per cent jLeinert et al.U l9931. 

In our survey we are sensitive to orbits with a period be- 
tween approximately 50 yr and 50,000 yr. For most of the sys- 
tems in our sample we find an orbital period of the order of a 
few thousand years. The literature data that we used are sev- 
eral decades old at most. The angular separation and position 
angle are not expected to have changed significantly over this 
time. If a star with significantly different values for p and/or <p 



is discovered with respect to the known companion star, this 
star is considered to be new. Several close binary systems are 
observed. The shortest orbital periods that we find (assum- 
ing circular and face-on orbits) are those for the companions 
of HIP62026 (57 yr), HIP62179 (68 yr), HIP76001 (84 yr), 
HIP64515 (87 yr), HIP80461 (106 yr), HIP62002 (142 yr), and 
HIP67260 (213 yr). Differences in angular separation p and 
position angle <p are therefore expected with respect to previ- 
ous measurements. For example, the observed star HIP76001 
has two companions with angular separations p\ and p2, and 
position angles <pi and <pz, respectively. In our AO survey we 
meas ure (pupi) = (0.2 5", 3.2°) and (p 2 ,n) = (1.48", 124.8°), 
while iTokovinir] i 19971) quotes (p u ipi) = (0.094", 6°) in 1993 
and (p 2 ,f2) = (1-54", 130°) in 1991. Taking the changes in p 
and <p into account, we determine whether an observed com- 
panion star was documented before or whether it is new. 

Table [3] summarizes the current status on binarity for 
Sco OB2. The number of presently known companion stars 
is listed as a function of primary spectral type and detection 
method. 

We define all companions that have only been detected 
with radial velocity studies as spectroscopic. All companions 
of which the presence of the companion has only been derived 
from astrometric studies of the primary are classified as astro- 
metric. All other companions, i.e. those that have been detected 
with (AO) imaging, Speckle techniques, interferometry, etc., 
are classified as visual. Note that our classification refers only 
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to available observations and not to intrinsic properties of the 
binary or multiple systems. 

All candidate companions that are found with our AO sur- 
vey are by definition visual. The overlap between the observed 
companions in our sample and the spectroscopic and astromet- 
ric binaries in literature is zero, although HIP78809 is flagged 
as "suspected non-single" in the Hipparcos catalogue. From 
our comparison with literature we find 41 new companions (14 
in US; 13 in UCL; 14 in UCL), while 33 of the candidate com- 
panion stars were already documented. 

Several of the stars currently known as 'single' could be 
unresolved systems. For a subset of the member stars, multi- 
ple companions have been detected using different techniques. 
Most systems with early spectral type primaries are found with 
spectroscopic methods. Amongst the A and F type stars, most 
Hipparcos member stars are single, although a fraction of these 
could be unresolved systems (see § 15.21 . 

Of the Sco OB2 Hipparcos member stars 37 are now known 
to be triple. Sixteen of these consist of a primary with two vi- 
sual companions. Fifteen consist of a primary, a visual com- 
panion, and a spectroscopic companion. These are only found 
in UCL and LCC. Six consist of a primary, a spectroscopic 
companion and an astrometric companion. Five quadruple sys- 
tems are known: two systems consisting of a primary and three 
visual companions (HIP691 13 and HIP81972) and two consist- 
ing of a primary, two visual companions and a spectroscopic 
companion (HIP80112 and HIP78384). HIP77820 is the only 
quintuple system, consisting of a primary, three visual compan- 
ions and one spectroscopic companion. The largest known sys- 
tem in Sco OB2 is HIP78374, which contains a primary, four 
visual companions and two spectroscopic companions. 

5.2. Binary statistics 

The binarity properties of a stellar populati on are usually quan - 
tified in 'binary fractions' (e.g., Reip urth & Zinneckerl [l993). 
Two common definitions in use are the multiple system fraction 
Fm and non-single star fraction F^s (since 1 - Fns is the frac- 
tion of stars that is single). Another frequently used quantity is 
the companion star fraction Fq 4 , which measures the average 
number of companion stars per primary star. These quantities 
are defined as 

F M = (B + T + ...)/ (S + B + T + .. .); (6) 
F NS = (2B + 3T + ...)/ (S +2B + 3T + .. .); (7) 
F c = (B + 2T + ...)/ (S + B + T + .. .), (8) 

where S , B, and T denote the number of single systems, binary 
systems and triple systems in the assocation. Table @] shows 
several properties of the three subgroups, including the binary 
statistics that are updated with our new findings. 

Figure^]shows the multiple system fraction Fm as a func- 
tion of spectral type of the primary for the three subgroups. For 
example, 86% of the B0-B3 type Hipparcos member stars of 
UCL have one or more known companion stars. Information 

4 Note that in lKouwenhoven et alJ l2004albl) we used an incorrect 
definition of F c ■ 



D Age S B T >3 F M F m F c 
(pc) (Myr) 



US 


145 


5- 


-6 


63 


46 


7 


3 


0.47 


0.67 


0.61 


UCL 


140 


15- 


-22 


131 


68 


17 


4 


0.40 


0.61 


0.52 


LCC 


118 


17- 


-23 


112 


54 


13 





0.37 


0.57 


0.45 


all 








303 


171 


37 


7 


0.41 


0.61 


0.52 



Table 4. Multiplicity among Hipparcos members of 
Sco OB2. The columns show the subgroup names (Upper 
Scorpius; Upper Centaurus Lup us; Lower Centaurus 
Crux), their distances (see de Zeeuw et al., 1999), the ages 
dde Geus. de Zeeuw & Lubl dl989t) : IPreibisch et all d2002l) for 
US: iMamaiek et all fcoda for UCL and LCC), the number 
of known single stars, binary stars, triple systems and N > 3 
systems, and the binary statistics (see § I5.2> . 



about the number of companion stars and the spectral type of 
the companion stars around the primary is not shown in the 
figure. The new companion stars found in our AO survey are 
indicated with the darker shaded parts of the bars. A trend be- 
tween multiplicity and the spectral type of the primary seems 
to be present, but this conclusion may well be premature when 
observational biases are not properly taken into account. Our 
detection of 41 new close companion stars shows that this is at 
least partially true: the number of A-type member stars which 
is in a binary/multiple system in US has doubled as a r esult of 
our su rvey. This strongly supports the statement made in Brownl 
(2001) that the low multiplicity for A- and F-type stars can at 
least partially be explained by observational biases. 

6. Conclusions and discussion 

We carried out a near-infrared adaptive optics search for com- 
panions around 199 (mainly) A and late-B stars in the nearby 
OB association Sco OB2. Our sample is a selecti on of the 
Hippa rcos membership list of Sco OB2 in Ide Zeeuw et alJ 
Jl999h . We find a total of 151 stellar components around the 
target stars. We use a simple brightness criterion to separate 
candidate companion stars (K$ < 12 mag) and probable back- 
ground stars (Ks > 12 mag). The validity of this criterion is 
verified in several ways (§|3}. Of the detected components, 77 
are likely background stars, and 74 are candidate companion 
stars. 

The 74 candidate companions occupy the full range in Ks 
magnitude, down to Kg = 12 mag. At the distance of Sco OB2, 
an M5 main sequence star has a Ks magnitude of approxi- 
mately 12 mag. The angular separation between primary and 
companion ranges from 0.22" to 12.4". These values corre- 
spond to orbital periods of several decades to several thousands 
of years. 

The J, H, and K s magnitudes for all observed objects are 
corrected for distance and interstellar extinction to find the ab- 
solute magnitudes My, M H , and M Ks . The subset of the com- 
ponents with multi-color observations are plotted in a color- 
magnitude diagram. All (except one) candidate companion 
stars are positioned close to the 5 Myr (for US) and 20 Myr 
(for UCL and LCC) isochrones. The background stars are po- 
sitioned far away from the isochrones. 
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For all observed primaries and companion stars the mass 
is derived from the absolute magnitude Mk s ■ The mass of the 
A and late-B primaries ranges between 1.4 M and 7.7 M Q and 
the mass of their secondaries between 0. 1 M Q and 3.0 M . Most 
observed companion stars are less massive than their primaries. 
Most of the systems with previously undocumented companion 
stars have a mass ratio smaller than 0.25. The mass ratio dis- 
tribution for the observed objects peaks around q — 0.1 and 
decreases for higher mass ratios. The minimum and maximum 
value for the mass ratio of the companion stars that we ob- 
served are q = 0.022 + 0.008 and 1.0 + 0.34. For systems with 
mass ratio q w 1 we cannot say with absolute certainty which is 
the most massive and therefore the primary star. The mass ratio 
distribution for binaries with late-B and A type primaries fol- 
lows the distribution f(q) = q~ r , with T - 0.33. This is sim ilar 
to the mass ratio distribution observed by S hatskv & Tokovinrr] 
(2002) for systems with B type primaries. Relatively few sys- 
tems with low mass ratio are found, excluding random pairing 
between primaries and companion stars. The uncertainty in the 
age of the UCL and LCC subgroups does not effect the mass 
ratio distribution significantly. 

A cross-check with visual, astrometric and spectroscopic 
binaries in literature shows that 41 of the candidate companion 
stars are new: 14 in US; 13 in UCL and 14 in LCC. The other 33 
candidate companions are already documented in literature. We 
analyze the presently known data on binarity and multiplicity 
in Sco OB2, including the 41 new companions that are found 
in our survey. We conclude that at least 41% of all Hipparcos 
member stars of Sco OB2 are either double or multiple, and 
we find a companion star fraction Fq = 0.52. These values are 
lower limits since the dataset is affected by observational and 
selection biases. 

Our AO observations are close to completeness in the angu- 
lar separation range 1" < p < 9". Next to the 199 target stars 
we find 50 companion stars in this angular separation range. 
This corresponds to an Fq of 0.25 in this range of angular sep- 
aration for A and late-B type stars in Sco OB2. This value (not 
corrected for incompleteness) is slightly higher than the val- 
u es for B type stars in Sco O B2 (Fq = 0.20 per decade of 
o: IShatskv & T okovinin. 2002) and pr e-main sequence st ars in 
Sco OB2 (F c = 0.21 per decade ofp: lKohler etalll2000h . 

At small angular separation (p < 3.75") no components 
other than the target stars are detected for magnitudes 12 mag < 
Ks < 14 mag. This 'gap' cannot be explained by observational 
biases or low-number statistics. A fully populated mass spec- 
trum for the companions (assuming random pairing of binary 
components from the same underlying IMF) is also incompati- 
ble with this gap. This implies a lower limit on the companion 
masses of ~ 0.1 M . This is consistent with our finding that the 
mass ratio distribution points to a deficit of low-mass compan- 
ions compared to the random pairing case. On the other hand, if 
we assume that the sources fainter than Kg « 14 mag are actu- 
ally physical companions, a gap may be present in the compan- 
ion mass distribution. We will carry out follow-up multi-color 
AO observations to further investigate this issue. 

The gap described above might indicate a brown dwarf 
desert, such as ob served for solar-type stars in the so- 
lar neighbourhood (Duauennov & Mavor, 1991). The pres- 



ence of this gap could support the embryo-ejection forma- 
tion scenario for brown dwarfs (e . g.lReipurth & Clarkeil200lt 
iBate. Bonnell. & Bromml 120031: iKroupa & Bouvierl 12003). 
This is further supported by the detection of 28 ca ndidate free- 
floatin g brown dwarf members of Sco OB2 by iMartfn et alJ 
ll2004l) . 

The observed Fm decreases wit h decreasing primary mass. 
Part of this trend was ascribed by Brown (2001) to observa- 
tional and selection biases. We lend support to this conclu- 
sion with our new AO observations. In particular, the multiplic- 
ity fraction for A-type Hipparcos members of US is doubled. 
Nevertheless our knowledge of the present day binary popula- 
tion in Sco OB2 is still rather fragmentary. None of the mul- 
tiplicity surveys of Sco OB2 so far has been complete due to 
practical and time constraints, but also because of a lack of full 
knowledge of the membership of the association. In addition, 
each observational technique used in these surveys has its own 
biases. For example, visual binaries are only detected if the an- 
gular separation between the components is large enough with 
respect to the luminosity contrast, which means that with this 
technique one cannot find the very short-period (spectroscopic) 
binaries. Moreover the observational biases depend on the way 
in which a particular survey was carried out (including how 
background stars were weeded out). 

We therefore intend to follow up this observational study by 
a careful investigation of the effect of selection biases on the in- 
terpretation of the results. This will be done through a detailed 
modeling of evolving OB associations using state-of-the-art N- 
body techniques coupled with a stellar and binary evolution 
code. The synthetic OB association will subsequently be 'ob- 
served' by simulating in detail the various binary surveys that 
have been carried out. This modeling includes simulated pho- 
tometry, adaptive optics imaging and Hipparcos data, as well as 
synthetic radial velocit y surveys. Examples o f this t ype of ap- 
proach can be found in Portegies Zwart et al. (2001) (for pho- 
tometric data) and Oui st & Linde gren (2000) (for Hipparcos 
data), where the latter study very clearly demonstrates that an 
understanding of the observational selection biases much en- 
hances the interpretation of binarity data. Finally, the time de- 
pendence of the modeling will also allow us to investigate to 
what extent stellar/binary evolution and stellar dynamical ef- 
fects have altered the binary population over the lifetime of the 
association. 

This combination of the data on the binary population in 
Sco OB2 and a comprehensive modeling of both the associ- 
ation and the observations should result in the most detailed 
description of the characteristics of the primordial binary pop- 
ulation to date. 
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Fig. 3. Angular separation between the target stars and ob- 
served components as a function of the component's K$ magni- 
tude (top) and component's Ks magnitude relative to that of the 
target star (bottom). The figures show the 41 new (dots) and 33 
previously known (plusses) companion stars (Ks < 12 mag), 
respectively. The triangles are the 77 background stars (Ks > 
12 mag; §[3}. The solid curves show the estimated detection 
limit as a function of angular separation. These figures clearly 
demonstrate that faint close companions are not detected be- 
cause of the extended PSF halo of the target star. The effec- 
tive field of view is 19.1" x 19.1". Companion star stars with 
angular separations larger than | V2 • 19.1" = 13.5" cannot 
be detected. For small angular separation (p < 3.5") no faint 
(Ks > 12 mag) components are found. This could indicate a 
gap in the companion mass function, or a lower mass limit for 
companion stars at small separation (see § 14.31 . 
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Fig. 4. The number of objects per angular separation bin ver- 
sus angular separation (bin size 2"). The distribution of pre- 
sumed background stars (dotted histogram) is clearly differ- 
ent from that of the companion stars (solid histogram). The 
expected angular separation distribution for background stars 
(in units of stars per bin size) is represented with the solid 
curve. The curve is normalized such that the number of stars 
corresponds to the expected number of background stars with 
12 mag < Ks < 16 mag (cf. Figure [3J. The excess of close 
(p < 6") companions could be caused by faint companion 
stars that are misclassified as background stars. Objects with 
p < 0.1" are too deep in the PSF halo of the primary to be de- 
tectable. Objects with p > 15" cannot be measured due to the 
limited field of view of the ADONIS/SHARPII+ system. 



Fig. 5. The expected background star distribution as a func- 
tion of angular separation and the Ks magnitude distribution of 
background stars as observed by Shatsky & Tokovinin (2002). 
The distribution of background stars over the field of view 
is assumed to be homogeneous. The triangles represent the 
observed stellar objects that we classify as background stars. 
The contours indicate the background star density in units 
of arcsec^'mag -1 . The distribution is normalized, so that for 
14 mag < Ks < 15 mag and 5" < p < 13" the number of 
background stars equals 19, i.e. the observed number of back- 
ground stars with these properties. Observational biases are not 
considered. The dashed line represents our criterion to sepa- 
rate companion stars (Ks < 12 mag) and background stars 
(Ks > 12 mag). The excess of background stars at small an- 
gular separations could be caused by faint companion stars that 
are misclassified as background stars. 
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Fig. 6. Top: the distribution of presumed companion stars over 
angular separation and position angle. Previously known com- 
panions and new companions are indicated with the plusses and 
dots, respectively. Bottom: the distribution of background stars 
over angular separation and position angle. The distribution of 
companion stars and background stars over position angle is 
random for angular separation smaller than 9.6" (see text). The 
companion stars are more centrally concentrated than the back- 
ground stars. 



Fig. 8. The distribution of the 74 detected candidate companion 
stars (solid histogram) and 77 presumed background stars (dot- 
ted histogram) as a function of Galactic latitude. The dashed 
histogram represents the distribution of Hipparcos member 
stars over Galactic latitude. The bin size is 5° for the three his- 
tograms. The distributions are clearly different: the companion 
star and Hipparcos member star distribution peak at about the 
central Sco OB2 latitude. As expected, the background star dis- 
tribution peaks at the location of the Galactic plane. 




■ . 



SilllS 




Fig. 9. Left: The hierarchical triple HIP68532. Both compan- 
ion stars (indicated with the white dots) were previously un- 
documented. Right: HIP691 13, another multiple system, where 
both close companion stars (white dots) were previously un- 
documented. A third known companion star with angular sepa- 
ration of 28" is not detected in our survey because of the large 
angular separation. 
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Fig. 11. Top: The relation between mass and Mk s magnitude. 
The 5 Myr isochrone (solid curve) is used to find masses of 
the companion stars in US, while the 20 Myr isochrone (dot- 
ted curve) is used for the UCL and LCC subgroups. The gray- 
shaded area is enclosed by the 15 Myr and 23 Myr isochrones, 
and represents the uncertainty in the age of the UCL and LCC 
subgroups. For objects at the mean Sco OB2 distance of 130 pc, 
Mk s = 6 corresponds to Ks = 11.6. Bottom: The derivative 
of the mass-magnitude relation. The conversion from Mk s to 
mass is most accurate where the absolute value of the deriva- 
tive dM/dM Ks is small. 



Fig. 12. Top: The mass distribution of the 74 companion stars 
(solid histogram) and corresponding primaries (dashed his- 
togram) observed in our AO survey. The bin size is 0.5 M . 
Masses are derived from their Ks magnitude and age (see 
§ 14.31 1. The median errors in the mass are 0.4 M Q for the pri- 
maries and 0.1 M for the companions. Bottom: Companion 
star mass versus angular separation. New and previously 
known companion stars are indicated as dots and plusses, re- 
spectively. Low-mass close companion stars are not detected 
due to the PSF wings of the corresponding primary star. 
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Fig. 13. Companion star mass versus primary star mass. The 
dashed lines represents the q = 0.25, q = 0.5, and q = 1 bi- 
naries and are shown to guide the eye. The 41 new companion 
stars and the 33 previously known companion stars are indi- 
cated with dots and plusses, respectively. The observed com- 
panion stars have masses lower than their associated primaries. 
Most systems with previously undocumented companion stars 
have q < 0.25. 
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Fig. 14. Top: The mass ratio distribution for the 74 systems for 
which we observe companion stars in our near-infrared AO 
survey (histogram). The mass ratio is defined as q — M c /M p 
where M c is the companion star mass and M p the mass of 
the corresponding primary star. All primaries are Hipparcos 
member stars with A and late-B spectral type. The q distri- 
bution does not change significantly if ages of 15 Myr or 
23 Myr for UCL and LCC are assumed (gray-shaded areas). 
The mass ratio distributions f(q) = q~ r are represented by the 
curves. For our observations we find F = 0.33 (dotted curve). 
IShatskv & TokovininN2002l) find F = 0.5 (dashed curve). The 
distribution which follows from random pairing (dash-dotted 
curve) is clearly excluded. The curves are normalized using the 
observed mass ratio distribution for the companion stars in the 
range 0.5 < q < 1.0. Bottom: The observed cumulative mass 
ratio distribution (histogram). The dotted and dashed curve are 
cumulative distribution functions corresponding to F = 0.33 
and F = 0.5, respectively. The fitted function described in 
Equation |4] is shown as the thin solid curve. The dash-dotted 
curve represents the cumulative q distribution for random pair- 
ing. 
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Fig. 15. The fraction of stellar systems which is multiple ver- 
sus the spectral type of the primary, for the three subgroups 
of Sco OB2. Only confirmed Hipparcos member primaries are 
considered here. The binarity dataset consists of literature data 
and our near-infrared AO survey. The light and dark grey parts 
of the bars correspond to literature data and the new data pre- 
sented in this article, respectively. The spectral types of the 
companion stars (not included in this plot) are always later 
than those of the primary stars. Apparently, the multiplicity is a 
function of spectral type, but this conclusion may well be pre- 
mature when observational biases are not properly taken into 
account. That this is at least partly true, is supported by our de- 
tection of 41 new close companion stars at later spectral types. 
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